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( 1 ) 


Personnel 


Dr. Li Kao Lee continues to be supported under this contract. 

(2) Expenditures 

The cumulative amount, $31,392.00, or 46% of the contract budget 
has been expended by the end of this the third quarter. 

(3) Technical Accomplishments 


3.1 Concentration Dependence of the Interdiffusion Coefficient. 

It now seems certain that the proper coefficient of diffusion to be 
used in the theory of Ostwald ripening of a two component system is the 
interdiffusion coefficient. This is the diffusion coefficient measured 
with respect to the center of volume of the system. This diffusion 
coefficient is conveniently determined with the diaphragm cell. 


In our first quarterly report, we showed that the concentration 
dependence of the interdiffusion coefficient can be easily measured using 
this cell if the volumes V, and V ? of the top and bottom compartments are 
eq_ual. In that report, we showed chat the interdiffusion coefficient, 
D(c), evaluated at c, the arithmetic mean of the top and bottom concentra- 
tions satisfied the equation 


in 


c ’ ( t ) - c" (t) " 
c '( 0 ) - c" (0) 


2AD(c)t 


(3.1.1) 


where i and A are the thickness and cross-sectional area, respectively, 
of the cell; V=V,=V~; and c'(t) and c"(t) are the bottom and top concen- 
trations, respectively, at time t, while c'(0) and c"(0) are the analagous 
concentrations at t=0. All other quantities being known, if c'(t) and 
c"(t) are followed as a function of t, D(c) can be determined directly 
from Eq ( 1 ) . 


We have found a new and more direct derivation of Eq (1), which exploits 
a Taylor series expansion of D(c) about c = c rather than c = 0 as done 
previously. The new derivation is much shorter and more direct. It leads 
to a new formula for estimating the systematic error, e, associated with 
Eq (1) (which is approximate). We find 

/ 1 \ 


D ( 2 ) (c)(ac ) 2 


D(C) 


x 100% 


(3.1.2) 


( 2 ) — 

where ac = c'(0) ^ c"(0) and '(c) is the second derivative of D(c) 
evaluated at c = c. Using reasonable values for the quantities in Eq (2), 
we find 


e - 1 .6% — (3.1.3) 

which is well within the random error ordinarily encountered when using 
the diaphragm cell. 

A research paper on this subject was delivered at the European Materials 
Research Society Meeting at Strasbourg, France May 13-15, 1985 and will be 
published in Journal de Physicue. A copy is attached. 



3.2 Grain Size Distribution in Ostwald Ripening as a Function of 
Volume Fraction of Grains 

We consider the precipitation of a single solute from a two component 
solution. Define the following quantities 

3 

c(“) - Equilibrium solubility in molecules/Cm of the solute in 
the solution. 

v - molecular volume of the solute in a precipitate grain in 
Cnr /molecule. 

D - Interdiffusion coefficient of solute and solvent in the 
parent phase in Cm /sec. 

y - Surface tension between a precipitate grain and the parent 
solution in erg /Cm. 

kT - Product of Boltzmann's constant times the absolute temperature, 
T, in ergs. 

R - Radius in Cm of a given precipitate grain. 

It is helpful to define some intermediate quantities. Let 

O = 2-yv/kT (3.2.1) 

Define a dimensionless time, t, 

t = tDvcH /a 2 (3.2.2) 

and a dimensionless radius, a, 

a = R/a (3.2.3) 

We find that the grain distribution, F(a,x) (where F(a,-r)da is the number 
of grains with radii between a and a + da) is given by 

F(a, t) = t' 4/3 F q (z) (3.2.4) 

Here 

z = a/x 1/3 (3.2.5) 


and 


FqU) = C Q z 2 (z o -z)' 0 (z+3/z 2 )' 6 exp(-6/(z 0 -z)), z<z Q 


where 


e = 2 + 


+ 61 
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(3.2.6) 

(3.2.7) 
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(3.2.9) 



■ The cut-off radius, z Q , satisfies the quartic equation 
ezj + 2z 3 - 3e 2 z 0 2 - 6ez 0 -3 = 0 (3.2.10) 

where 

e = (4.55 9 ) 1/2 (3.2.11) 

and cp is the precipitate grain volume fraction which occurs once the 
steady state of the Ostwald ripening has been established. The normaliza- 
tion constant, C Q , depends upon <p. 

Figure 1 shows the function F (z) for various values of <p. These curves 
were computed as follows: A value°of <p was chosen and e was calculated 

using Eq (3.2.11). Eq (3.2.10) was solved for z using this value of e. 

Next, the exponents e, 6, and 6 were computed using Eqs (3.2.7) - (3.2.9). 

The value of the normalization constant was determined from the integral 

z Q 

jdz z 3 F Q (z) = 1 (3.2.12) 

o 


which guarantees the conservation of solute mass. 

Figure 1 demonstrates that competition for solute between solute 
grains has a substantial effect on F (z), even for small values of <p. 
The grain size distribution becomes broader and the maximum 3 radius 
(cutoff radius) increases with <p. Because the area under z'T (z) is 
conserved (See Eq (3.2.12)), the brodened distributions for cp°> 0 are 
also shorter. Detailed derivations of Eqs (3.2.1) - (3.2.12) will be 
presented in the final report. 
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Abstract - We derive an approximate analytical technique for extracting the 
concentration dependence of the interdiffusion coefficient from diffusion 
measurements carried out with the diaphragm cell. The systematic error 
incurred in using the technique is estimated and found to be no greater 
than the random error ordinarily encounted in these experiments. 


I - INTRODUCTION 

In microgravity, Ostwald ripening in phase separating immiscible liquids can be 
observed without interference from convection and sedimentation. Two component 
systems such as succinonitri le/water or succinonitri le/ethanol are easily studied, 
because they are optically transparent. In order to compare results obtained with 
these' systems to the theory of Ostwald ripening /1.2.3/, the interdiffusion coef- 
ficient, and the interfacial energy and densities of the separating phases must be 
measured. 

If the two components form an ideal solution, the interdiffusion coefficient is 
independent of composition and can be measured directly in earth's gravity using 
the diaphragm cell /4/. In the case of the two monotectic systems, NC (Chp) ? CN/H ? 0 
and NClCH^lpCN/C^H^On, mentioned above, however, the interdiffusion coef T 1 cient is 
a funetiofi of composition. In as much as the diaphragm cell yields the integral of 
the interdiffusion coefficient over a range of concentration, the integral data must 
be unfolded to obtain the desired quantity /5/. 

Below, we show that under certain conditions, the concentration dependence of the 
interdiffusion coefficient can be obtained directly from the data without unfolding. 
There is a systematic error incurred in this procedure, but it is of second order 
and no larger than the random! error ordinarily incurred with the cell. 

II - THEORY OF THE MEASUREMENT 

The diaphragm cell consists of two volumes, V, and V ? , separated by a sintered glass 
diaphragm disk of cross sectional a r ea, A, an6 nominal thickness, t /4./. The sen- 
sitivity of cell depends upon the ratio, A ft. Because of the torturous nature of 
the diffusion paths among the glass particles in the diaphragm, i cannot be directly 
observed. Rather, the cell constant, kit. is determined by observing the diffusion 



coefficient of a solution (often aqueous KC1) whose interdiffusion coefficient is 
absolutely known. 

This being done, we assume that the cell is loaded at t=0 with a solution of the two 
components of concentration c'(0) occupying V. and concentration c"(0) occupying 
V~. (To prevent convection, the cell is ordinarily oriented so that i is parallel 
t6 the gravity vector, and the less dense solution is above the more dense.) Thus 
take c * ( 0) > c"(0). If the volumes V. and V ? are well stirred, a concentration 
gradient, dc/dx, exists only across the diaphragm. The gradient is related to the 
flux, J, through Pick's first law 

J = - D ( c ) dc (1) 

dx 


where the interdiffusion coefficient Ibl is allowed to be a function of the concen- 
tration, c. Eq. (1) can also be written, Jdx =-D(c)dc, and integrated on the left 
side from 0 to £ and on the right side from c’ to c" to obtain 

c" 

Je = - J D(c)dc (2) 

c 1 


By conservation of mass, the rate of loss of solute from V., -V^dc'/dt, must exactly 
equal the rate of gain of solute, Vpdc"/dt by Vp. Hence, 


V~dc" 

Z 3t 


JA 


-V.dc ' 
1 3t 


(3) 


Eqs. (2) and (3) can be combined to obtain 


d(c ' -c" ) 

cTE 



( 4 ) 


As written, Eq. (4) cannot be integrated directly without knowledge of the func- 
tional form of D(c), which is ordinarily not known. 

In what follows, we shall show that Eq. (4) can be integrated accurately after- a 
Taylor series development of D(c), provided that V 1 = V^V- In tins circumstance, 
conservation of solute mass dictate_s that the mean concentration c = (c 1 4 c")/2, 
is independent of time and if c' = c 4 ac, then c" = c - ac. Define the variable 
y = c - c and transform the integral, 1, on the right hand side of Eq. (4) to 

c 1 AC 

I = J D ( c ) dc = I D("c 4 y) dy ORIGINAL PA(sE*JS (5) 

c" -ac OF. POOR QUALfTY 


Expand D(c 4 y) about D(c) 

D(c 4 y) = I 
n=0 

where D n ( c” ) = (d n D(c )/dc n ) 
term by term. The result 

1 ; ! D r^> 

n=0 TnTTT ! 


in a Taylor series 
n ! 

|c = c. Substitute 
i S 


((AC) 


ne i 



in y, 


Eq. (6) 


) 


into Eq. 


( 6 ) 


(5) and integrate 


(7) 


When written to third order in Ac, Eq. (7) reads 

, ORIGINAL PAGE* IS 

I = 2D 0 (c)ac ♦ (l/3)D 2 (c)(Acr OF POOR QUALITY (8) 
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Note that the term of order (Ac) cancels out. Factoring Eq. (8) and using 
Ac = ( 1/2) (c'-c") , we have 


I - D(c) (c ' -c" ) 

1 + (l 

\ D 2 (c ) ( AC ) 2- | 


. V 

’ ) D(cf) J 


(9) 


where the "zeroth^ derivative D 0 (c) = D(c). The integral, 1, thus can be accurately 
represented by D(c)(c'-c") assuming that the systematic error, c. 


1 \ D 2 (c)(ac) 2 

*)'~i 


is small. 

We replace 1 by D(F) (c 1 -c")and integrate Eq. (4) to obtain 


£n 


( c'(t) - c"(t)\ _ 
y c ' ( 0 ) - c ,r TUT ) = 


2AD(c) t 
fV 


(ID - 


Eq. (11) is identical with the usual result /4 , 7/ , except V = V. = V~ and the 
interdiffusion coefficient is a function of the mean concentration, L c. If the 
concentrations, c'(t) and c"(t), are monitored as a function of time, then a plot 
of the left_hand side of Eq. (U) as a function of t, yields a straight line with 
slope -2AD(c )/ fV from which 0(c) can be extracted if V and A ft are known. 


III - DISCUSSION 

If_Ac is not too large, we can safely assume that D~Cc)(ac) is a small fraction of 
D(c), say 1/10. In this circumstance, using Eq.(10j,.we estimate c = 1/60, which 
is the same order of magnitude as the random error ordinarily incurred with the 
diaphragm cell. To test Eqs. (10) and (11), a range of solutions of increasing ac 
should be prepared and tried with the cell. So long as for each solution D(c) 
determined using Eq. (11) is independent of ac, it can be concluded that the sys- 
tematic error is small. In as much as the diffusion will be slow when ac is smell, 
special care should be taken to keep the cell isothermal over the course of the 
experiment. 
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